Comment on "Effects of spatial dispersion on electromagnetic 
surface modes and on modes associated with a gap between two 

half spaces" 

G. L. Klimchitskaya^ and V. M. Mostepanenko^ 
^North-West Technical University, Millionnaya St. 5, St. Petersburg, 191065, Russia 
"^Noncommercial Partnership "Scientific Instruments", 
Tverskaya St. 11, Moscow, 103905, Russia 

Abstract 

Recently Bo E. Sernelius [Phys. Rev. B 71, 235114 (2005)] investigated the effects of spatial 
dispersion on the thermal Casimir force between two metal half spaces. He claims that incorporat- 
ing spatial dispersion results in a negligible contribution from the transverse electric mode at zero 
frequency as compared to the transverse magnetic mode. We demonstrate that this conclusion is 
not reliable because, when applied to the Casimir effect, the approximate description of spatial 
dispersion used is unjustified. 

PACS numbers: 71.10.-w, 73.25.+i, 71.36. +c, 68.47.De 
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Ref. [l| is devoted to the question of how spatial dispersion affects the Casimir force 
between real metal plates at nonzero temperature. As is mentioned in Ref. there are two 
quite different predictions in the literature regarding the thermal Casimir force between real 
metals. According to one prediction at separations of a few hundred nanometers 

the thermal correction to the Casimir force is negligible and agrees in a qualitative sense 
with the ideal metal result. At large separations, where the magnitude of the total force 
is equal to the thermal correction (z > Q nm at T = 300 K), the results for real and ideal 

n n 

metals coincide. According to another prediction [51, |6|, at short separations below one 
micrometer the thermal correction for real metals can be as large as 10-20% of the total 
Casimir force. At large separations the thermal force between real metals is one-half of the 
ideal-metal result. Mathematically the distinction between the two proposed formalisms 
lies in the fact that in Refs. there is a nonzero contribution from the transverse 

electric mode at zero frequency, whereas in Refs. this contribution is equal to zero. A 
nonzero result for the contribution from the transverse electric mode at zero frequency is 
obtained by using the dielectric function of the plasma model {2, 3] or the Leontovich surface 
impedance Q]. A zero_result follows by application of the ^rude dielectric function wi_th a 



relaxation parameter 



,l6|. As was demonstrated in Refs. j?, 8], the prediction of Refs. js,^] 
is excluded by experiment at 99% confidence within a wide separation region from 300 to 
500 nm, whereas the predictions of Refs. are consistent with experiment. There is 

also the qualitative disagreement between the prediction of Refs. and the results of 

the first modern measurement of the Casimir force (9| at z ~ 1 um. At 2; = 1 um the large 
thermal correction, predicted in Refs. [5|, 16| , achieves 19% of the force in ideal metal case, 
whereas the experimental uncertainty in jsj was 5-10%. The predicted thermal effect was 
not. however, observed in g Ref. Q elaims that the melusion of spatial dispersion leads 
to a result which is almost identical to that obtained from the Drude model in Refs. [5|, l6|. 
Below we demonstrate that including spatial dispersion into the theory of the Casimir force 
as in Ref. is unjustified and, thus, the result obtained in this reference cannot be trusted. 

To find the electromagnetic surface modes associated with an empty gap between the 
two half spaces, Ref. ij] assumes in Sees. Ill and IV the validity of the standard continuity 
boundary conditions 
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where n is the normal to the boundary directed inside medium, and the subscripts n, t refer 
to the normal and tangential components, respectively. Here E is the electric field strength, 
D is the electric displacement, B is the magnetic induction, and only nonmagnetic materials 
are considered so that for the magnetic field H = B (the subscript 1 refers to a vacuum 

n 

and subscript 2 to a semispace). To describe the spatial dispersion, Ref. [1[ assumes that 
the longitudinal and transverse dielectric functions depend on both the wave vector and 
frequency eti = Ski{q,uj). However, in the theory of the thermal Casimir force taking into 
account spatial dispersion, both assumptions used in Ref. i.e., the boundary conditions 
([T]) and the dielectric permittivities depending on both q and u, are unjustified. 

While on the subject of boundary conditions, we recall that in Casimir-type problems 
the Maxwell equations of an electromagnetic field in a metal can be written as 
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rotE + — ^ = 0, divD = 0, 

c ot 
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rots ^ = 0, divB = 0. 
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Here it is assumed, that there are no current and charge densities of the external sources 
(i.e., given independently of E, D and B). The contribution from the conduction electrons 
is taken into account by the relation for the determination of D 

dD BE ^ . 

where i is the volume current density induced by the fields E and B. 

In electrodynamics with spatial dispersion the physical fields E and B at the interface 
are usually finite, though the displacement D can tend to infinity [10]. Taking this fact 
into account, and integrating the Maxwell equations ([2]) over the thickness of the boundary 



layer (see, e.g., Ref. lUj), one reproduces the first two conditions of Eq. ([T]), but arrives to 



modified third and fourth conditions 
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Elt = E2t, Bin = B2nj (4) 

D2n-Din = Ana, [nx {B2- Bi)] = —j. 

c 



Here the induced surface charge and current densities are given by 

"2 1 r2 



If 1 /" dD 

a = — / div[n x \D x n\]dL j = — / —:—dl. (5) 



In linear electrodynamics without spatial dispersion the material equation connecting D 
and E takes the form 

Dk{r,t) = I eu{r,t-t')Ei{r,t')dt' (6) 

J — oo 

(here we consider a medium with time- independent properties). According to Eq. ([2]), the 
electric displacement D{r,t) at a position r and time t is determined by the electric field E 
at the same point r (the spatial dispersion is absent) but at different times t' <t (generally 
speaking, the temporal dispersion is present). It can be easily shown that the substitution 
of Eq. ([6]) in Eq. ([5]) leads to a = 0, j= 0. In other words when only a temporal dispersion 
is present the boundary conditions (jll) coincide with the standard continuity conditions 
([T]). However, it is unjustified to use the boundary conditions ([T]) in the presence of spatial 
dispersion as is done in Ref . Ill . Specific examples illustrate that in this case neither a nor j 

□ n 

is equal to zero [10|, |13[. Note that the boundary conditions (j4]) are the exact consequence of 
the Maxwell equations with sharp boundary surfaces and they are obtained for macroscopic 
physical fields. It is not correct to mix them up with the boundary conditions arising in 
perturbative theories and for fictitious fields (see below). ^ 

We turn now to the second assumption used in Ref. pj], i.e., to the description of spa- 
tial dispersion in terms of dielectric functions depending on q and u. If only a temporal 
dispersion is present, one can substitute the Fourier transformations 

/oo /"OO 
E{r,uj)e~''^'duj, D{r,t) = D{r,uj)e-''^'duj (7) 
OO J — oo 

in Eq. ([6]) and arrive at 

Dk{r, uj) = Ekiir, uj)Ei{r, u), (8) 

where 

poo 

eki{r,u)= / ikiir,T)e'''^dT (9) 

is the frequency-dependent dielectric permittivity and r = t — t'. Eq. and bound ary 
conditions ([T]) are used in all traditional derivations of the Lifshitz formula 
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I20I describing the van der Waals and Casimir force between parallel plates (half 
spaces) possessing only a temporal (in other words "frequency" ) dispersion. 

If the material of the plates possesses both temporal and spatial dispersion, Eq. IQ should 
be generalized to 

Dk{r,t) = j dt' j dr'iki{r,r',t-t')Ei{r',t'). (10) 



If the medium were uniform in space (i.e., all points were equivalent), the kernel e would 
not depend on r and r' separately, but only on the difference R = r — r' . In this case, by 
performing the Fourier transformation 

/oo /» 
du / dqE{q,uj)e'^'''-^'^ (11) 
-oo J 

(and the same for D), and substituting it in Eq. (fTOj) . one could introduce the dielectric 
permittivities 

/oo /• 
dr J dRhi{R^T)e-''^^'^-'^^\ (12) 

These permittivities depend on both the wave vector and frequency and bring Eq. ffTOl) into 
the form 

Dk{q, uj) = Ekiiq, uj)Ei{q, u) (13) 

in analogy to Eq. 

In Ref . pj] , however, the system under consideration is not uniform due to the presence of a 
macroscopic gap between the two half spaces. Because of this, it is impossible to assume that 
the kernel e depends only on R and r, and hence it is not possible to introduce £ki{q, ^^)- In 
fact in the presence of boundaries the kernel i for systems with spati^ dispersion depends 
not only on R and r but also on the distance from the boundary [10]. An approximate 
phenomenological approach to dealing with this situation is described in Ref. [l^. For 
electromagnetic waves with a wavelength A the kernel e{r, r', r) is significantly large only in 
a certain vicinity of the point r with characteristic dimensions a ^ A (in fact for nonmetallic 
condensed media a is of the order of the lattice constant). One can then assume that e is a 
function of R=r r', except for a layer of thickness a adjacent to the boundary surface. If 
one is not interested in this subsurface layer, the quantity eki{q,u>) may be used to describe 
the remainder of the medium. 

This approximate phenomenological approach is widely applied in the theory of the 
anomalous skin effect (see, e.g., Ref. 2l|) for the investigation of bulk physical phenom- 



ena involving electromagnetic fields. In Ref. 



2l| some kind of fictitious infinite system is 



introduced and the electromagnetic fields in this system are discontinuous on the surface. 
This discontinuity should not be confused with the discontinuity of physical fields of a real 
system in the presence of spatial dispersion given by Eqs. (jlj) and ([5]). (There is, however, 
another approach which describes the anomalous skin effect in polycrystals in terms of the 
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local Leontovich impedance taking into account the shape of Fermi surface 22|]). The fre- 
quency and wave vector dependent dielectric permittivity in the presence of boundaries is 
also successfully applied in some other problems, e.g., in the theory of radiative heat transfer 



231] or in the study of electromagnetic interactions of molecules with metal surfaces 24 1. 



All these applications are in more or less good agreement with experiment. The boundary 
effects can be taken into account by the boundary conditions (jlj) supplemented in some way 
by so called "additional boundary conditions" . 

It is unlikely, however, that the phenomenological approach of Ref. [21I using the fre- 
quency and wave vector dependent dielectric permittivity would be applicable for the cal- 
culation of Casimir force between metal surfaces, where the boundary effects for vacuum 
electromagnetic oscillations are of prime importance and contribute significantly to the to- 
tal result. Ref. [l| uses the quantities Ekiiq^uj) in a spatially nonuniform system of two half 
spaces separated by a gap without justification, and determines the surface modes from 
the continuity conditions ([T]) instead of the boundary conditions (jlj) valid in the presence 
of spatial di=pe..io„. As was emphasized in Ret. M, in the case ot a bounded medium it 
is inadmissible to use the same kernel in Eq. ( |TOl) which is appropriate for an unbounded 
medium. According to Ref. [25| , this leads to a violation of the law of conservation of energy. 
It is then not surprising that in doing so Ref. arrives to an incorrect contribution from 
the transverse electric mode to the Casimir free energy. However, in addition to Ref. |l| 
some other recent publications (see, e.g., Refs. [26|, l27|]) apply the dielectric permittivities 
eki{q,u!) obtained for an unbounded medium to the theory of the Casimir effect which es- 
sentially depends on the presence of boundaries. 

One more shortcoming of Ref. jl| is that it uses the usual Lifshitz formula for the free 



energy derived in the presence of only temporal dispersion Q, Q, Q, 0, Q, 
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20|, to 



compute the effects of spatial dispersion. It has been known, however, that with the inclusion 
of spatial dispersion the free energy of a fluctuating field takes the form = J^l + AJ^, where 
J-'l is given by the Lifshitz formula derived in a spatially local case and written in terms 
of the Fresnel reflection coefficients, and AJ-" is an additional term which can be expressed 
in terms of the thermal Green's function of electromagnetic field and polarization operator 



ii 

28l | . The review paper [281] contains a few references to incorrect results by different authors 



obtained by the substitution of Ski{q, 00), taking account of spatial dispersion, into the usual 
Lifshitz formula. During the last years several more similar papers have been published 
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(see, e.g., Refs. |26l. l27j). It can be true that the Lifshitz formula written in terms of general 
reflection coefficients is applicable in both spatially local and nonlocal situations. However, 
as far as the exact reflection coefficients in a spatially nonlocal case are unknown, the use 
of some approximate models, elaborated in the literature for applications different than the 
Casimir effect, may lead to incorrect results for AjF and create serious inconsistencies with 
experiment. 

To conclude, the main purpose of Ref. [ll "to find out how spatial dispersion affects the 
Casimir force between real metal plates" is in our opinion not achieved. Instead of using 
the boundary conditions (jll) valid in the presence of spatial dispersion, and the appropriate 
generalization of the Lifshitz formula, Ref. |l| uses the continuity conditions ([T]) and usual 
Lifshitz formula derived in the presence of only temporal dispersion. Ref. pj] also applies the 
phenomenological dielectric permittivity which depends on frequency and wave vector to the 
calculation of the Casimir force, a phenomenon primarily determined by virtual photons. 
These unjustified assumptions are responsible for the main conclusion of Ref. ll that spatial 
dispersion results in practically the same "dramatic effect" that was previously obtained in 
Ref. 51] from the Drude dielectric function: the contribution from the transverse electric 
mode at zero frequency "completely vanishes if dissipation is included", and is negligibly 
small compared to the contribution from the transverse magnetic mode if dissipation is 
absent. Ref. [l| does not inform the reader that theoretical approaches with zero or negligibly 
small contribution from the transverse electric mode at zero frequency have already been 

n n 

excluded experimentally at 99% confidence in Refs. [7|, 18|. By now such approaches are 
rejected by 4 different experiments 29|] including the experiment of Ref. 91]. At the same 
time, the Leontovich impedance approach, which leads to a nonzero contribution from the 
transverse electric mode at zero frequency, is both consistent with experiment and takes 
dissipation into account , |30| . The physical reasons why the Leontovich impedance 



approach is to be preferred are articulated in Ref. [30]. 

As to the role of spatial dispersion in the Casimir effect, it was first investigated using the 
Leontovich surface impedance in Ref. [311 • The Leontovich approach can be considered as 
an approximate phenomenological approach. As was shown above, the approach of Ref 
in spite of the author's claims to the contrary, is also phenomenological and approximative. 
At the moment there is no fundamental theory which describes the thermal Casimir force 
between real metals incorporating spatial nonlocality. Because of this, one should choose 
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between different approximate aproadies, and in so doing tfie role of experiment must not be 
underestimated. It is well understood that at laboratory temperatures of about 300 K and 
at Casimir plate separations of a few hundred nanometers (the frequency region of infrared 
optics) spatial dispersion does not play any role. It can be confidently neglected also in the 
region of the normal skin effect, i.e., at separations between the plates greater than 4-5 //m. 
If some approximate phenomenological approach leads to the opposite conclusion, there are 
strong reasons to beheve that it is not rehable. 
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